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The hnear confinement in quarkonium is generahsed as a minimal tree, 
with interesting geometrical properties. This model binds tetraquarks more 
easily than the additive model used in earlier investigations. 
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1. Introduction 

The stability of compound systems is a delicate issue. In atomic physics, 
binding H~ cannot be established using any Hartree type of wave func- 
tion /(ri)/(r2). The positronium molecule, Ps2 = (e+, e"*", e~, e~) is rather 
weakly bound below the two-atom threshold. The configuration (p,e~^,e~) 
is unstable, but it is stabilised by a second electron, leading to the positron- 
ium hydride, below the H + Ps threshold. There are many other examples. 

In molecular physics, the dimer ^He2 is barely bound, while the higher 
^IIe„ clusters with n > 3 have much deeper binding per atom, illustrating 
the Thomas effect [IJ. In the same He-He potential, the lighter ^He2 dimer 
is unbound, even for spin singlet. It was attempted to calculate ^He„ with 
n = 3, ... 10, which were found unbound. Every reasonable physicist would 
have given up, but fortunatel}{3 some colleagues were curious enough to push 
the calculation to higher n, and discovered that ^He„ is bound for n > 35, 
or less if there are also some ^He admixed [2]. 

In nuclear physics, some improbable configurations turn out eventually 
stable, due to the pairing interaction. For instance, while ^He is unstable. 
^He is stable. Seen as a (a, n, n) three-body system, it is Borromean, as 
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none of the two-body subsystem is bound. Even more remarkable is that, 
though "^He is unstable, one finds again stability for ^He, something like a 
small drop of neutron star around two protons. 

In hadron physics, too, exotics probably require a sophisticated arrange- 
ment of the quarks and antiquarks. After the wave of baryonium candi- 
dates, which were not confirmed, and the absence of firmly-established light 
pentaquarks, some scepticism now prevails about multiquark spectroscopy. 
However, many sectors remain to be explored, for instance, mesons with 
charm 2 or with beauty —2. 

The best known category of multiquarks includes the deuterium and 
other nuclei, but, due to the repulsive core of nuclear forces, each nucleon 
keeps it identity. The Yukawa mechanism is by no means restricted to the 
nucleons, and acts between several other hadron pairs, provided they contain 
light quarks. In particular, the pion-exchange is attractive in the DD* + 
D*D channel. As compared to the deuteron, the potential is weaker but 
experienced by heavier constituents. Hence this approach was considered as 
successful when the X(3872) was discovered just above the DD* threshold. 
However, there is no hard core here, and one cannot escape the short-range 
forces, which presumably result from a direct interaction between the heavy 
and light quarks. For a recent review, and refs. to the pioneering papers, 
see, e.g., g]. 

Some years ago, an interesting coherence has been discovered in the 
colour-spin operator entering the spin-spin interaction known as chromo- 
magnetism Some expectations values are larger (and with the right sign 
for inducing attraction) for some multiquark configurations than for the sum 
of the decay products. This is the case in particular for the l:i{uuddss) as 
compared to A-l-A, or for the 1987-vintage pentaquark [5j, (Qqqqq), as com- 
pared to (Qq) + (qqqq), where (qqqq) denotes (uuds) or permutations. The 
first estimates were carried out with unbroken flavour SU(3)f symmetry, and 
assuming that the quark-quark correlation is the same for multiquarks as 
for ordinary hadrons. Unfortunately, breaking the flavour symmetry ben- 
efits more to the threshold than to the H or P and thus spoils bindings. 
Also multiquarks are more dilute systems than ordinary hadrons. Thus the 
strength of spin-spin forces is appreciably reduced. Altogether, the H and 
the P seem likely unbound when all effects are taken into account. 

If the chromomagnetic interaction, once rather promising, turns out dis- 
appointing, it is natural to address the possibility of binding using the chro- 
moelectric interaction, in particular through its properties under symmetry 
breaking and its A^-body character (this is not a sum of pairwise terms). 
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2. Flavour independence and symmetry breaking 



Symmetry breaking is known to lower the ground state of any HamiL 
tonian. Tlie simplest example is + + \x in one dimension, with the 
lowest energy ^(A) = 1 - AV4 < -E(O). More generally if 



H = Hn + XHi 



(1) 



where Hq is even under some symmetry, and Hi odd, then the variational 
principle with the even ground state solution of the A = case as a trial 
wave function gives 



E{X) < E{0) . 



(2) 



However, in most cases, the threshold energy benefits more from this effect, 
and symmetry breaking deteriorates stability. This is what happens with 
the H vs. A + A in the chromomagnetic model, when SU(3)f is broken. 
Another example is the breaking of permutation symmetry in molecules. 
The equal mass case, say (^■^, ^u^, //~), a rescaled version of the positro- 
nium molecule, is bound by about 3% below the threshold for dissociation 
into two neutral atoms (internal annihilation is neglected here). Consider 
now the asymmetric configuration m+, M~ ,m~), with fi~^ average of 
and m~^. The Hamiltonian can be written as 



H{M- 



, m 



, m 



) = H{n~^,n'^,iJ, ) + 



1 

4M 



1 

4m 



-pI+pI- 



-pD 



(3) 

demonstrating that E(M~^m~^ , M ,m ) < E{iJ,'^ , fi^ , fj, , fi ). However, 
meanwhile, the threshold evolves from 2£'2(/U,/^) to E2{M,M) + E2{m,m) 
which is lower. Indeed, it has been shown from accurate four-body calcula- 
tion that stability holds only for 1/2.2 < M/m < 2.2. See, e.g., [6]. 

However, a miracle occurs when breaking charge conjugation. Consider 
a mass configuration (M, M, m, m) and a potential that does not depend 
on the masses, as the Coulomb potential in atomic physics, or the static 
potential in QCD (fiavour independence). Then 



H{M, M, m, m) = H{^, /x, ^, /i) + 



1 

4M 



1 

4m 



(pI+pI 



pI 



■pI), (4) 



Again, the ground-state energy is lowered by the odd term, but the thresh- 
old energy remains constant, 2 E{M,m) = 2 E2{fJ,, fJ-), as the two-body en- 
ergy depends only on the reduced mass. Thus stability is improved. For 
molecules, it is, indeed, observed that the relative excess of energy as com- 
pared to the threshold, evolves from about 3% for M = m to nearly 17% 
for M ^> m, e.g., for the hydrogen molecule. 
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Similarly, explicit quark model calculations with a flavour-independent 
potential indicate that stability is reached for {Q,Q,q,q) against dissocia- 
tion into two flavoured mesons if the mass ratio is large enough. 

Another miracle is that while the spectroscopy of exotic hadrons is usu- 
ally a matter for hot controversy, the stability of {Q,Q,q,q) in the large 
M/m limit has reached an overall consensus. See, e.g., O |8] for refs. This 
is also confirmed in a QCD sum-rule calculation [9]. 

The question now is whether the double charm is heavy enough to get 
stability. On the experimental side, there are candidates for the double 
charm baryons from the SELEX experiment, but they are not confirmed in 
other experiments, yet, and double-charm mesons have never been searched 
for. See, e.g.,[10j and refs. there. On the phenomenological side, explicit 
four-body calculations have been carried out of (Q, Q, q, q) and similar con- 
figurations, using potential models tuned to reproduce ordinary mesons and 
baryons. The usual conclusion is that double charm is not sufficient to bind, 
and that double beauty would be safer. However, Rosina et al. |7J, pushing 
very far a variational calculation with a realistic potential, got the state 
(c, c, c, d) weakly bound. 

3. Steiner tree of confinement 

The problem, however, is whether the interquark potential can be safely 
extrapolated from mesons and baryons to multiquarks. The usual prescrip- 
tion is the colour-additive rule 

y(ri,r2,...) = -^^K~\jv{rij) , (5) 

i<j 

where v{r) is the quarkonium potential. This is perhaps justified for short- 
range contributions corresponding to colour-octet exchange. But there is 
not a serious reason, that it should be adequate for confinement. In fact, 
many years ago, the generalisation of v{r) = r (in units where the string 
constant is unity) has been proposed by several authors. For refs., see, e.g., 
[8] . It has been recently supported by detailed simulations using the lattice 

QCD [n]. 

For baryons, the extension is the so-called y-shape potential 

V3{vi,V2,vs) = min((ii + d2 + ds) , (6) 

s 

where di is the distance of the i^^ quark located at Vi {i = 1, 2, 3) to a 
junction s whose location is adjusted to minimise V3. This corresponds to 
the well-known problem of Fermat and Torricelli to link three points with 
a minimal network. See Fig. [TJ 
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For tetraquarks, the potential, also pictured in Fig. [U reads (with dij = 

\\ViVj\\) 

U = min{di3 + d2i, du + ^23, V4} , 

F4 = min ( ||uiSi|| + ||f2Sl|| + IIS1S2II + ||s2f3|| + ||s2f4|| ) , 

Sl,S2 

with the minimum of the two possible quark-antiquark links, sometimes 
referred to as the "flip-flop" model, and the connected flux tube, itself 
minimised by varying the location of the Steiner points si and S2- 




Fig. 1 . Generalisation of the linear quark-antiquark potential of mesons to baryons 
(left) and to tetraquarks, where the minimum is taken of the flip-flop (centre) and 
Steiner tree (right) configurations. 

A first investigation using the potential ([7]) concluded to the absence of 
exotics [12]. However, Vijande et al. [13] used a more systematic variational 
expansion of the wave function and in their numerical solution of the four- 
body problem, found a stable tetraquark ground state. Moreover, unlike 
[12\ , they considered the possibility of unequal masses, and found that sta- 
bility improves if the quarks are heavier (or lighter) than the antiquarks, in 
agreement with the earlier argument about favourable symmetry breaking. 

Recently, an inequality has been derived for this tetraquark potential, 

^<^(IN| + ||y||) + ||z|| . (8) 

The four-body Hamiltonian with this upper bound splits into three simple 
one- variable Hamiltonians with linear confinement, and can be solved an- 
alytically. In the limit of large M/m, the stability below the dissociation 
threshold is recovered [8]. 

4. Outlook 

The recent progress in understanding confinement has inspired more 
realistic quark potentials which have been used for studying multiquark 
spectroscopy. Remarkably, this new modelling of confinement gives bet- 
ter attraction than the conventional colour-additive models, and thus pre- 
dicts a richer spectrum of exotics. This is confirmed by a detailed study of 
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(Q, Q, g, q) configurations. It is our intend to apply the same dynamics to 
pentaquark and dibaryon configurations. 
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